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1. INTRODUCTION 

Numerical method is widely used in computer science. It helps researchers, engineers, and students 
to approximate the solution of complicated problems so that the final solution consists of only addition, 
subtraction and multiplication operations. Numerical method forms an important part of solving Initial Value 
Problems (IVPs) in Ordinary Differential Equations (ODEs), especially in cases where the exact solution is 
difficult to obtain orin the absence of closed form analytic formula. 

The Euler method is a first-order numerical procedure for solving ODEs with a given initial value. 
It is an effective [1] and the easiest method [2] to solve IVPs in ODEs. Euler method has been use in variety 
field such as in communication and security to get better solution. Regarding [3], Euler method has better 
performance in modeling uncertainties. It is suppeoted by [4] where as experimental results indicate that it is 
much simpler and faster, more sensitive and robust. Although Euler provides simple procedure for solving 
ODEs, there are issues that have compelled the use of other, more complex, methods such as Runge Kutte. 
One of them is “accuracy” as reported by [5-6]. Improvements to the Euler method have attracted much 
attention including [1-2], [5-7], resulting in modified Euler methods. 

This paper is aim to discover a new algorithm as accurate as possible with exact solution by using 
average concept. The new algorithm proposed in this paper name it as Cube Arithmetic. This idea to develop 
this new algorithm is by enhance other algorithm proposed by Qureshi [7] and Nurhafizah [8]. Qureshi used 
average of arithmetic mean for two points at coordinate y=[f(xn,yn)+f(xn+lyn+1/2. The equation is 
referring as Euler Arithmetic shows improvement in accuracy and speed compared to Euler Method. [8] uses 
concept of Qureshi [7] but choose Harmonic mean and Contraharmonic mean in the equation. 
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Thus, this paper proposes a modified Euler method using the Cube mean concept. The method was 
constructed by extracting the Euler Arithmetic (EA) method in [9] and using Cube mean in the equation. 
Propose algorithm is compared between modified Euler known as Euler Arithmetic with the exact solution. 
The accuracy of the proposed method is demonstrated by the comparing error results of with step sizes of h 
0.1, 0.01 and 0.001. 


2. PROPOSED ALGORITHM 

The algorithm is based on the Cube Mean in Euler Arithmetic. Euler Arithmetic (EA) is a well- 
known technique for improving the Euler method. Euler methods used in [1] and [2] were extracted as basis 
for the algorithm. Cube mean within two coordinate points of function were used to improve the methods. 
Equation (1) is a basic formula of the Euler method. 


Yast = Yn tAtf (fo, Yo) (1) 


Equation (2) by [1] was modified using the average concept. 
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The proposed average in cube mean for the two midpoints is written as Equation (3) 
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Modifications to the slope of the function at estimated midpoints of (xo. Vo) and (5 y,) would 


improve the stability and accuracy of Euler. 


3. RESEARCH METHOD 
There are a few processes to present a mathematical model. First, it can be presented by 
mathematical software and algorithm development. Before using the mathematical software, algorithm 
development must be done first. Algorithm refers as instructions for solving problem using simple 
language [10-11]. According to [12], construction process in mathematical software includes as follows. 
a) The design and analysis algorithms 
b) Algorithm coding 
c) Testing 
d) Details documentation 
e) Distribution and maintenance of the software 
Algorithm that has been developed will be converting to the programming language. In this study, 
the code program is written in SCILAB 6.0 to test the effectiveness of the algorithm. The effectiveness of the 
algorithm will be identifying by taking itaccuracy and time of taken to process the algorithm. This paper will 
focus on accuracy by comparing the exact solution and solution between proposed method known as Cube 
Euler Arithmetic and modified Euler naming as Arithmetic Euler. 
Cube Arithmetic (CA), the proposed algorithm in Equation (3), was represented using SCILAB 6.0. 
SCILAB is a powerful, low-cost [13]-[14] and flexible tool for researchers, engineers, and students to use in 
mathematical computer application.Scilab was developed by the French National Institute for Informatics 
and Automation Research (INRIA) [15].The effectiveness of CA was determined through the comparison of 
both CA and modified EA solutions, focusing on accuracy and processing time. 


4. RESULTS AND ANALYSIS 

Results of three first order ODEs were recorded. These results were obtained by testing CA with 
three sets of linear ODEs using three different step sizes: 0.1, 0.01 and 0.001. These sets were tested for 
efficiency using the minimum and maximum error for the whole cycle. Table | illustrates the exact solutions 
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for the ODEs. Table 2 illustrates the comparison results between EA and CA methods. Relative error based 
on [17] was calculated as: 


Error= [Ex Ev] , Ex=Exact_value and Ev=Euler’s_modified_value 
Table 1. Set of Problem First Order ODEs 
Equation _ Exact Solution Initial Values Interval of Integration Source 
y’=-0.5y en y(0)=0 0<x<20 [16] 
y=y ew y(0)=1 0<x<20 [16] 
y’=-10y el) y(0)=1 0<x<20 [16] 
Table 2. Results of Errors Using Various Step Size h 
Method Euler Arithmetic Cube Arithmetic 
Step Size h=0.001 h=0.01 h=0.1 h=0.001 h=0.01 h=0.1 
Problem | 0 0 0.000083 0 0.000001 0.000076 
Problem 2 0 0.000003 0.000355 0 0.000003 0.00302 
Problem 3 0 0.000401 0.117900 0 0.000300 0.029020 


The results for each method were compared to the exact solution using maximum error. It is evident 
from Table 2 that the proposed CA method provided a more accurate result than EA at higher h sizes. 
Smaller step size produced similar results for both methods. 

Problem | demonstrated that the results of the CA method is better than EA at A=0.lwith a 
maximum error of 0.000076 and 0.000083 respectively. Figure | illustrates the differences between the two 
methods. Results at h=0.01 showed that CA scored 0.000001 compared to EA at 0. Both methods scored the 
same at h=0.001. Problem 2 demonstrated that the results for both CA and EA methods are almost the same 
at h=0.1, 0.01 and 0.001, as shown in Figure 2. 
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Figure 1. Differences between the two methods 
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Figure 1. The results for both CA and EA methods 
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Problem 3 demonstrated that the CA method scored better results compared to EA at each step size. 
At h=0.1, CA scored 0.029020 whilst EA scored 0.117900, and at h=0.01, CA scored 0.000300 whilst EA 
scored 0.00040, as shown in Figure 3. 
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Figure 2. The results for both CA and EA methods 


It can be summarized thatthe proposed CA method provides a more accurate result than EA in 
solving ODEs at higher step size. 


5. CONCLUSION 

This paper has proposed Cube Arithmetic, a new method using modified Euler as finding of this 
study. The feasibility of the proposed algorithm was compared with Euler Arithmetic and tested using 
SCILAB. Subsequently, the Cube Anthmetic compared with Euler Arithmetic scheme with the exact 
solution. Usually, the ordinary Euler method using a small step size gives the solution almost to the exact 
solution. However, in this study has been demonstrated that Cube Arithmetic provided solutions that are 
similar to exact solutions at small step size and also at higher step size. The benefits of using small step size 
will give a higher accuracy. Therefore a higher step size would reduce complexity and processing time. As a 
conclusion, Cube Arithmetic can be used as an alternative algorithm to solve ODE problems. 
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